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ABSTRACT 

We present Keck/DEIMOS spectroscopy of stars in the recently discovered Milky Way satellites 
Hydra H, Pisces H, and Laevens 1. We measured a velocity dispersion of km s“^ for Pisces H, 

but we did not resolve the velocity dispersions of Hydra H or Laevens I. We marginally resolved the 
metallicity dispersions of Hydra H and Pisces H but not Laevens 1. Furthermore, Hydra H and Pisces H 
obey the luminosity-metallicity relation for Milky Way dwarf galaxies (([Fe/H]) = —2.02 ± 0.08 and 
—2.45±0.07, respectively), whereas Laevens 1 does not (([Fe/H]) = —1.68±0.05). The kinematic and 
chemical properties suggest that Hydra H and Pisces H are dwarf galaxies, and Laevens I is a globular 
cluster. We determined that two of the previously observed blue stars near the center of Laevens 1 
are not members of the cluster. A third blue star has ambiguous membership. Hydra H has a radial 
velocity (vheiio) = 303.1 ± 1.4 km s“^, similar to the leading arm of the Magellanic stream. The 
mass-to-light ratio for Pisces 11 is 370^240 ^q/Lq- It is not among the most dark matter-dominated 
dwarf galaxies, but it is still worthy of inclusion in the search for gamma rays from dark matter 
self-annihilation. 

Subject headings: galaxies: dwarf — Local Group — galaxies: abundances 


1. INTRODUCTION 

Ultra-faint dwarf galaxies (UFDs) harbor a wealth 
of information about dark matter and nucleosynthesis. 
They are the most dar k matter-dominated objects known 
(jSimon fc Gehal 120071 1. Their large dark-to-luminous 
mass ratios and their small sizes imply large central den¬ 
sities of dark matter. Hence, they are excellent tar¬ 
gets for the dete ction of gamma rays fr om dark matter 
self-annihilation (iBonnivard et al.ll2015ll . They also con¬ 
tain the larg est concentrations of metal-poor stars of any 
galaxy type ()Kirbv et al.ll2008bl : iFrebel et al.ll2014ll . The 
most metal-poor stars in UFDs could be the direct de¬ 
scendants of the first generation of stars in the Universe. 
If so, then their compositio ns are direct samp les of Pop¬ 
ulation HI nucleosynthesis ([Frebel et al.ll20I0ll. _ 

Th e Sloan Digital Sky Survey ISDSS. lAbazaiian et al.l 
1200911 revolutionized the field of dwarf galaxies a decade 
ago by discovering the first UFDs. In total, more 
than a dozen new UFDs were found in the SPSS (e.g., 
Wil^an et al.l[2005l : IZucker et al.ir2006l: iBelokurov et a,lJ 

200711 . In the span of roughly five years, SPSS and its 
successor, the Sloan Extensio n for Galactic Unde rstand- 
ing and Evolution (SEGUE, lYannv et al.ll2009l l. more 
than doubled the number of known Milky Way satel¬ 
lites. The field is experiencing another rejuvenation with 
the arrival of several new imaging surveys with deeper 
photometry and/or coverage of the Southern sky: the 
Dark Energy Survey (DES), the Survey of the MAgel- 
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lanic Stellar History ISMASH. iNidever et al.l[201^ . the 
Panoramic Sur vey Telescope fc Ra pid Response System 
(Pan-STARRS, IKaiser et al.l[20pl. A TLAS at the VLT 
Survey Telescope (iShanks et al.ll2015ll . and independent 
imaging with the Dark Energy Camera (DECam) at the 
CTIO/Blanco telescope. To date, DECam imaging, in¬ 
cluding DES, has enabled the discovery of 12 new Milky 
Way satellites (iKoposov et al.ll20I5a:lBecht ol et al.ll20I5l : 
iKim fc JerienI l2015allbl : iKim et al. _1201^ |bll . Pan- 
STAR RS discovered two satellites ([Laevens et al.l I20l4 
1201511 , and ATLAS and SMASH have each discovered one 
([Belokurov et al.l [20141: iMartin et aDl20I5r FI. Together, 
these surveys discovered 15 additional satellites in less 
than a year and a half. Their success can be attributed 
to deeper photometry and surveying a different part of 
the sky than SDSS. 

The photometric discovery is the first step in identi¬ 
fying a new dwarf galaxy. In order to classify a stel¬ 
lar system as a galaxy, it should exhibit some evidence 
for dark matter in the form of a lar ge mass-to-light ra¬ 
tio, m etallicity dispersion, or both (IWillman fc Strad^ 
1201211 . These criteria can be tested only with spec¬ 
troscopy. Two of the newly discovered satellites. Retic¬ 
ulum 11 and Horologium 1, have been spectroscopically 
confirmed as dwa rf galaxies by both velocity and metal¬ 
licity dispersions (iSirnon et al.ll201^ iWalker et al.ll201^ 
iKoposov et al.l[2015bll . The other systems have tenta¬ 
tive classifications as galaxies or globular clusters (GGs) 


based on their luminosities and half-light r adii (r/, ). GC s 
have Th < 30 pc regardless of luminosity (iHarrisl l IMi), 
but the Th of galaxie s increases with luminosity (e.g., 
IBelokurov et al.l l2007f) . The half-light radii of dwarf 
galaxies overlap with GGs {vh ~ 30 pc) at My ~ —2, 
and they grow to ru > 100 pc at My < —5. 


^ ILaevens et al.l II2014I) . using Pan-STARRS, discovered 
Laevens 1 simultaneously with IBelokurov et al.l fl2Q14) . who used 

ATLAS. 
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Thi s study conce r ns th ree recently discovered satel¬ 
lites. iMartin et al.l (|2015f) discovered Hydra II in DE- 
Cam images taken as part of SMASH. Its large half- 
light radius, 68 ± 11 pc, strongly suggests that it is a 
dwarf galaxy. Hydra H is especiall y interesting for its 
proximity to the Magellanic Clouds. IMartin et al.l raised 
the possibility that it is associated with the Large Mag¬ 
ellanic Cloud (LMC), which is po tentially true for many 
of the DES-discovered satellites (TKoposov et ahll^lSaJ: 

iBechtol et al.ll201^lDeason et al.l[2015l l. iBelokurov et al l 

(|20ir)H discovered Pisces H in SDSS images, and they 
confirmed an overdensity of main sequence turn-off stars 
with deeper KPNO/MOSAIC images. Its size, rt ~ 
60 pc, suggests that it is also a galaxy. However, 
no spectroscopy has b een ob tained since it s disco very. 
IBelokurov et alJ ()2ni4[ l and iLaevens et al.l (j2014fl co¬ 
discovered Laevens 1/Cratei0 with ATLAS and Pan- 
STARRS, respectively. Its half-light radius is only about 
20 pc, and its luminosity is My ~ —5. Therefore, it 
can be tentatively classified as a GC, but spectroscopy is 
required for a definitive identification. 

We observed Hydra H, Pisces H, and Laevens 1 with 
Keck/DEIMOS in order to identify them as galaxies or 
clusters. We completed the identification by measuring 
both velocity and metallicity dispersions. We also quan¬ 
tified the metallicities of a few stars in these galaxies to 
get a hint of their capability to enhance their own metal¬ 
licities. 

2. PHOTOMETRY AND ASTROMETRY 

In order to design DEIMOS slitmasks, we needed a 
catalog of coordinates for stars that were potentially 
members of the three satellites. We also needed colors 
and magnitudes o f stars to identify candidate members. 
ISand et HI (120121) published a photometric and astro- 
metric catalog for Pisces H, but no such published cat¬ 
alog exists for Hydra H or Laevens 1. We downloaded 
publicly available images for Hydra H, and we obtained 
new Keck/LRIS images of Laevens 1. 

2.1. Hydra II 

IMartin et al.l (I2015D discovered Hydra H in 
SMASH/DECam images. The Hydra H images, taken 
on 2013 March 20, are publicly available through the 
NO AO Science Archive^. We downloaded the calibrated 
images in the DECam g and r filters. These images are 
flat fielded, and they have astrometry headers. We used 
SExtractor (jBertin fc Arnouti 119961) to identify stars 
and measure their magnitudes in each of the g and r 
images. We discarded objects with class_star < 0.8 
in order to weed out galaxies and image artifacts. We 
corrected m agnitudes for extinct ion according to the 
dust maps of iSchlegel et al.l ()1998l SFD98). 

We selected stars for spectroscopy based on their col¬ 
ors and magnitudes. We drew a region in the color- 
magnitude diagram (CMD) in the approximate shape of 
the red giant branch (RGB) and horizontal branch (HB). 
Figure [T] shows this region in gray. We assigned priori¬ 
ties for spectroscopic selection within the region based on 
magnitude but not on color. These priorities were used to 


resolve conflicts where slitmask design constraints forced 
a choice among two or more stars. Brighter stars were 
given higher priorities. The stars that were able to be 
placed on the DEIMOS slitmask are shown as blue points 
(members) and red crosses (non-members). Figure [2] 
shows the sky coordinates of the spectroscopic targets. 

Our spectroscopic selection includes two potential HB 
stars, 189086 and 191385, and one asymptotic giant 
branch (AGB) star, 194563. Our kinematic membership 
selection (Section [5T|) includes all three stars. However, 
star 191385 is 0.4 mag above the HB. Therefore, we ruled 
it a non-member. Star 189086 has a color and magnitude 
consistent with the HB. Although it is the most distant 
member in our Hydra H sample, we kept it in our list of 
members. These decisions do not affect our results in a 
measurable way (Section [Oil . 

2.2. Pisces II 

ISand et Ml (j2012[) observed Pisces H with Magel¬ 
lan/Megacam in the g and r filters. We used their 
published astrometric and photometric catalog. They 
corrected for extinction using the SFD98 dust map. 
We drew a spectroscopic selection region around the 
RGB (Figure [T]). Because Pisces H has a smaller half¬ 


light radius (1(1, 

Sand et ?LL||2012tl than Hvdra II (1'7. 

IMartin et al.ll2015 

), the possible targets are denser on the 


sky, which causes more conflicts for spectroscopic selec¬ 
tion. As a result, the selection region did not include the 
HB so that RGB stars could be selected instead. (RGB 
spectra lend themselves more easily to the measurements 
of radial velocity and metallicity than HB spectra.) As 
for Hydra H, spectroscopic priority was given to brighter 
stars. 

2.3. Laevens 1 

D. Perley kindly observe d Laevens 1 for us with 
Keck/LRIS (jOke et al.lll995|) on 2015 March 23. Simul¬ 
taneous exposures were obtained for 130 s in the g filter 
in the blue arm and 120 s in the / filter in the red arm. 

The images were reduced with LP ip4!l, which provides 
flat fielding, astrometric solutions, and photometric cal¬ 
ibration. 

We identified s tars and measured their magnitudes 
with DAOPHOT (IStetsoTill1987l[2?frn) . The LRIS point 
spread function (PSF) was asymmetric, and it varied 
over the field. We allowed DAOPHOT to choose the 
best analytic PSF, which was a Penny function (a two- 
dimensional Gaussian core with Lorentzian wings). Ad¬ 
ditionally, there was a look-up table that allowed the PSF 
to vary linearly over the field. Again, we corrected for 
extinction using the SFD98 dust map. 

B ecause the half-light radius of Laevens 1 (0(5- 
0(6, IBelokurov et al.l [2014 ILaevens et al.l [20T3) is even 
smaller than Pisces H, we again selected stars on the 
RGB, not the HB. However, we also included three 
bright, blue stars near the center of t he system. The 
nature of these stars is controversial. IBelokurov et al.l 
(|2014D speculated that these are core helium-burning 
blue loop stars, which would imply that Laevens 1 has 
formed stars wi thin the last few hun dred Myr, making it 
a dwarf galaxy. ILaevens et al.l (|2014[ ) also identified these 

^ http://www.astro.caltech.edu/$\sim$dperley/programs/lpipe.html 


® We call the object Laevens 1, as is the convention for GCs. 
® http://www.portal-nvo.noao.edu/ 
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Figure 1. Color-magnitude diagrams of the three satellites. Blue points are spectroscopically confirmed members, whereas red crosses are 
non-members. Small black points are stars within 4' of the satellite’s center that we did not observe spe ctroscopical ly. The gray regions 
circumscribe the stars that were considered possible members. The ridgeline for the metal-poor GC M92 IIClemll20Q6l) is shown in cyan. 
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Figure 2. The sky coordinates of spectros copic targets. The values shown are displacements from the satellite center, as measured by 
[Martin et al.l II2015I1 . iBelokurov et al.l II2010I1 . and lLaevens et 33 II2014I1 . Blue points are spectroscopically confirmed members, whereas red 
crosses are non-members. Small black points are stars with go < 22 not targeted for spectroscopy. The DEIMOS slitmask outline is shown 
in green. Because the field for Laevens 1 is smaller than the other two satellites, the axis ranges are smaller. 


stars, but they did not favor their interpretation as blue 
loop stars because young, blue stars at fainter magni¬ 
tudes were absent. We obtained spectra of all three stars 
in order to determ ine if they are membe rs of Laevens 1. 
We also note that iBonifacio et al.l (|2015ll found a popu¬ 
lation of stars brighter and bluer than the main sequence 
turn-off but too faint for spectroscopy. If they are not 
blue stragglers, these stars could signify the presence of 
a Gyr old population in Laevens 1. 

For the purposes of measuring metallicities, it is 
convenient to have photometry in a uniform system. 
Therefore, we converted Cousins / magnitudes into 
SPSS i magnitud es following the conversion formula of 
I.Tordi et al.l ([200l): i = I + 0.21(E - I) + 0.34. The for¬ 
mula depends weakly on i? — / color, which we approxi¬ 
mated as 0.5 mag for all RGB stars. Due to the impreci¬ 


sion of this assumption, we added 0.05 in quadrature to 
the i photometric errors. This error floor corresponds to 
an error in i? — / color of 0.2 mag. For comparison, the 
full range of i? — / color for an old, metal-poor RGB is 
about 0.4 mag. 

3. SPECTROSCOPY 
3.1. Observations 

We obs erved one slitma sk for each satellite with 
DEIMOS (jFaber et al.l[200l on 2015 May 18. Table ffl 
lists the exposure times of each slitmask. We used the 
I200G grating with a ruling of 1200 lines mm“^ and 
a blaze wavelength of 7760 A. The grating was tilted 
such that the center of the CCD mosaic corresponded to 
7800 A. Slits were 0.7" wide. This configuration gives 
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Table 1 

DEIMOS Observations 


System 

UT Date 

targets Airmass 

Seeing 

Individual Exposures 

Total Exposure Time 





(") 

(s) 

(s) 


Milky Way Satellites 

M22 

2009 Oct 13 

64 

1.5 

0.6 

2 X 900 

1800 


2009 Oct 14 

64 

1.5 

0.6 

2 X 900 -1- 420 

2220 

Hydra II 

2015 May 18 

49 

1.6 

0.8 

4 X 1200 -1- 1 X 900 

5700 

Pisces II 

2015 May 18 

20 

1.8 

0.8 

3 X 1380 

4140 

Laevens 1 

2015 May 18 

14 

1.4 

0.9 

3 X 1380 

4140 

Laevens 1 (star 1717) 

2015 Mar 22 

1 

1.2 

0.6 

1 X 1800 

1800 

Radial Velocity Standard Stars 

HD 38230 

2012 Apr 19 


1.4 

1.0 

1 X 120 

120 

HD 151288 

2013 Apr 13 


1.1 

1.0 

1 X 80 

80 

HD 103095 

2013 May 5 


1.1 

0.7 

1 X 60 

60 

BD-18° 5550 

2013 May 18 


1.3 

0.8 

1 X 350 

350 

HD 88609 

2015 May 18 


1.2 

0.8 

1 X 240 

240 

HD 122563 

2015 May 18 


1.4 

0.8 

1 X 75 

75 

HD 187111 

2015 May 18 


1.2 

0.9 

1 X 80 

80 

BD+23° 3912 

2015 May 18 


1.0 

0.9 

1 X 300 

300 

HD 109995 

2015 May 19 


1.1 

0.9 

1 X 300 

300 

Telluric Standard Stars 

HR 4829 

2015 May 19 


1.1 

0.9 

1 X 60 

60 

HR 7346 

2015 May 19 


1.1 

0.7 

1 X 120 

120 


an approximate wavelength range of 6300-9100 A at a 
resolution of 1.3 A FWHM {R - 6000 at 7800 A). The 
exact wavelength range of each spectrum depends on the 
location of the slit on the slitmask. The starting and end¬ 
ing wavelengths of the spectra vary by ^ 300 A across 
the slitmask. We also obtained internal flat field and arc 
lamp exposures in the afternoon and morning for cali¬ 
bration. 

As mentioned above, iBelokurov et al.l (j2014D identified 
several bright, blue stars in the vicinity of Laevens 1. 
These stars would be unusual in an old GC or dwarf 
galaxy. Our slitmask included two of these stars. We 
observed another blue star, 1717, with a single 0.7" slit. 

We also observed nine radial velocity standard stars 
and two telluric (hot) stars with a long slit. The slit 
width was 0.7". We used a long slit that spanned the 
entire length of the DEIMOS field of view. This allowed 
us to refine the wavelength solution based on night sky 
lines over the entire CCD mosaic. Section 13.21 describes 
our approach to the wavelength solution. 

3.2. Reductions 

As in our previous papers (e.g.. iSimon fc Gehal 120071 : 
iSimon et al.ll2nill : iKirby et al.ll2013al 1201,^ . we reduced 
the DEIMOS spectra using a slightly modified version of 
the spec2d ID L data reduction pipeline developed by th e 
DEEP2 team (|CooDer et al.l[201^ iNewman et al.ll2013f) . 
We introduced a few updates to those procedures for this 
data set. The updates include improvements to the wave¬ 
length solution as determined from sky emission lines and 
corrections for the effects of differential atmospheric re¬ 
fraction. In the first stage of the sky line wavelength 
tweaking, a zero point offset is determined during the 
main reduction pipeline. We improved the tracing of the 
sky lines across each slit at low S/N, which is particu¬ 
larly important for our long-slit observations of bright 


template stars with short integration times. We also 
added a second stage to the sky line fitting: after the 
reductions are completed, a quadratic fit is performed 
to the sky line wavelengths in each extracted spectrum, 
and then the variation of each of the quadratic fit pa¬ 
rameters as a function of position on the slit mask is fit 
with a polynomial. This process corrects for errors in 
the flexure compensation system as well as temperature 
changes between the times at which the arc frames and 
the science frames were obtained (Geha et ah, in prep.). 

Atmospheric refraction shifts the position of the star 
light within the slit during the observations. Shifts in 
the spatial direction (along the slit) result in curvature 
of the object profile on the detector, and we changed the 
DEEP2 extraction algorithm to account for this effect 
(|Kirbv et al.ll201^ . Shifts in the wavelength direction 
(across the slit) result in velocity offsets that vary as a 
function of wavelength. Using the airmass at the time 
of observation, the angle between the slit and the par¬ 
allactic angle, and the measured seeing, we computed 
a correction for this velocity shift and applied it to the 
extracted spectra. 

In order to take maximum advantage of these improve¬ 
ments, we constructed a new empirical library of tem¬ 
plate stars for radial velocity measurements to replace 
the template set m ost DEIMOS dwarf g alaxy studies 
have employed since ISimon fc Geh^ (l2007f ) . We selected 
a sample of metal-poor stars spanning a range of effec¬ 
tive temperature, surface gravity, and metallicity, and 
we observed them by orienting the slit north-south and 
slowly driving the telescope such that the star moved 
steadily across the slit. This process ensured that the 
star light uniformly illuminated the slit during the expo¬ 
sure. To make sure that the template star observations 
contained strong enough sky lines for accurate adjust¬ 
ments to the wavelength solution (as described above), 
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Figure 3. Example spectra for one red giant in each of the satel¬ 
lites. The top two spectra have high S/N, and the bottom spectrum 
has low S/N . We fit model spectra (red) to the high-S/N spectra 
fSection 14.211 . Pisces II 10694 is carbon-rich with CN absorption 
visible between 8300 A and 8400 A. The model spectrum has a 
normal carbon abundance, so there is excess CN absorption in the 
observed spectrum compared to the model. 


we integrated for at least 60 s during each template obser¬ 
vation even if the star crossed the slit in a much shorter 
amo unt of time. We also replaced the telluric template 
from iSimon fc Gehal (1200711 with higher S/N observations 
using the same techniques. 

Figure |3] shows example I-D spectra of one member 
star in each satellite. The spectra shown for Hydra H and 
Pisces H have high S/N, while the spectrum for Laevens I 
has low S/N. Also shown are the best-fit model spectra 
(described in Section 14.21) for the stars with S/N high 
enough to measure metallicity. Because the S/N is low 
for the star in Laevens 1, we did not measure a metallicity 
for it, and no model spectrum is shown. 

4. SPECTROSCOPIC MEASUREMENTS 
4.1. Radial Velocity Measurements 

We measured radial velocities of stars in the three 
Milky Way satellites by comparing their spectra to the 
radial velocity template stars ’ spectra. In a manner sim¬ 
ilar to lSimon fc Gehal (1200711 . we optimized the DEEP2 
surve y’s redshift measurement technique (iNewman et al.l 
1201,411 for measuring stellar velocities rather than galaxy 
redshifts. This technique involves computing = 
(target flux — template flux)^/(flux error)^ for each tem¬ 
plate star. The relative velocity between the target 
and the template is shifted to find the minimum y^. 
The velocity of the star (uobs) is the one that mini¬ 
mizes y^ among all of the templates. As INewman et al.l 
pointed out, this technique is a ge neralization of a cross¬ 
correlation (|Tonrv fc DavisI [197911 that allows for differ¬ 
ent values of flux error for each pixel. 

Small astrometric errors or imprecision in the align¬ 
ment of the slitmask cause each star to fall at some dis¬ 


placement from the exact center of the slit. Mis-centering 
tra nslates to a shi f t in t he wavelength scale. Follow¬ 
ing [S^o^&(^eh3 (|2007ll . we corrected all of the stel¬ 
lar velocities to a standard geocentric frame by finding 
the velocity offset (ugeo) required to align the observed 
telluric absorption features with the features in the hot 
star spectra. The observed wavelength regions included 
in the velocity correction were 6866-6912 A (B band), 
7167-7320 A, 7593-7690 A (A band), and 8110-8320 A. 
The velocity corrections were determined with the y^ 
technique described in the previous paragraph. The ve¬ 
locities quoted in this paper are in the heliocentric frame: 
Vheiio = I'obs + Vgeo + Vcorr, where Ucorr is the conversion 
from the geocentric frame to the heliocentric frame. 

We estimated random and systematic uncertainties 
separately. To estimate the random uncertainty, we re¬ 
sampled the spectra for 10^ Monte Carlo trials. In one 
trial, the flux of each pixel was drawn from a Gaussian 
probability distribution. The mean of the distribution 
was the original flux value, and its variance was the flux 
variance estimated by spec2d. We re-measured Uheiio (in¬ 
cluding both Uobs and Ugeo) for all of the trials, and we 
took the random uncertainty, drand'C) to be the standard 
deviation of the Uheiio measurements. 

There is some minimum velocity uncertainty, even 
for noiseless spectra. This uncertainty can arise due 
to uncorrected flexure in DEIMOS, spectral template 
mismatches, errors in the wavelength solution, and un¬ 
known causes. We estimated SsysV by calculating the 
difference in Uheiio for two independent measurements of 
the same set of stars. The ideal sample for estimating 
6sysV will have many high-S/N spectra. Instead of using 
the dwarf galaxy sample, we used spectra of 52 red gi¬ 
ants in the GC M22. The spectra were taken with the 
same slitmask on consecutive nights (see Table [I]). The 
data from each night were reduced separately accord¬ 
ing to Section 13.21 We measured Uheiio and estimated 
Aandr' as described above. We then calculated the quan¬ 
tity (Ubelio,! 'Chelio, 2 )/\/Aand,l'C^ 4“ Aand, 2 r:^ T 2(5sys'C^ 
for each pair of observations of the same star. The 
standard deviation of this quantity should be 1 for well- 
determined errors. The value of SsysV required to satisfy 
that conditio n is 1.49 km s~^ . The final error bar for 
each star is \/(^rand'C^ + Ays'd^ • 

Our estimate of d^ysV ignores some potentially impor¬ 
tant sources of error. First, the repeat measurements 
were obtained with the same slitmask. As a result, some 
of the slit mis-centering—for example, due to astrometric 
errors in the slitmask design—would be repeated in the 
same way across the two consecutive nights. Second, the 
slitmasks were observed at about the same hour angle 
on each night. Although we applied wavelength correc¬ 
tions for flexure and differential atmospheric refraction, 
these errors would have been of the same sign and similar 
magnitude for both nights. Any residual, uncorrected er¬ 
ror would be similar across both nights. Third, the M22 
spectra have fairly high S/N. Consequently, the same 
template spectrum was assigned to 47 out of 52 M22 
stars. This fraction would be lower in spectra with lower 
S/N. 

All of these effects would potentially lower the mea¬ 
surement of SsysV artificially. In order to quantify the 
effect of a spuriously low value of dsys'c, we also calcu- 
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Table 2 

Target List 


ID 

RA (J2000) Dec (J2000) go {g — r)o (g — I)o 

S/N^^ 

'^helio 

Member? 

logs 

[Fe/H] 


(mag) (mag) (mag) 

(A-i) 

(km s“^) 

(K) 

(cm s“^) 



Hydra II 


196246 

12 21 27.50 

-31 57 22.6 

19.68 

0.63 


54 

161.8 ± 1.6 

N 




184647 

12 21 28.13 

-32 06 36.7 

19.75 

0.63 


46 

116.3 ± 1.9 

N 




199452 

12 21 30.18 

-31 55 58.0 

19.94 

0.85 


59 

5.1 ± 1.6 

N 




189524 

12 21 32.70 

-32 03 51.5 

21.90 

0.60 


8 

272.9 ± 7.9 

N 




188838 

12 21 33.19 

-32 04 29.4 

20.66 

0.67 


33 

319.0 ± 2.3 

N 




194103 

12 21 34.74 

-31 59 21.8 

20.93 

0.57 


25 

300.0 ± 4.2 

Y 

5095 

2.05 

-2.40 ±0.32 

193538 

12 21 35.93 

-31 59 54.0 

19.79 

0.67 


66 

304.2 ± 1.7 

Y 

4918 

1.46 

-1.95 ±0.12 

198021 

12 21 36.31 

-31 56 38.1 

21.40 

0.68 


19 

63.9 ± 3.4 

N 




195726 

12 21 36.74 

-31 57 54.1 

21.13 

0.60 


22 

306.5 ± 3.5 

Y 

5036 

2.11 

-2.48 ±0.33 

190646 

12 21 38.18 

-32 02 38.2 

20.23 

0.61 


44 

86.6 ± 1.8 

N 




194563 

12 21 38.51 

-31 58 56.4 

20.57 

0.48 


34 

302.9 ± 3.3 

Y 




191521 

12 21 39.68 

-32 01 48.8 

19.91 

0.96 


45 

202.5 ± 1.8 

N 




191385‘> 

12 21 40.08 

-32 01 57.7 

20.70 

0.16 


9 

316.9 ± 19.9 

N 




192206 

12 21 40.37 

-32 01 12.1 

19.97 

1.02 


68 

28.0 ± 1.5 

N 




193286 

12 21 40.96 

-32 00 07.0 

21.27 

0.54 


14 

299.9 ± 5.4 

Y 




194920 

12 21 41.02 

-31 58 34.4 

21.58 

0.54 


16 

294.5 ± 5.9 

Y 




194405 

12 21 41.05 

-31 59 04.1 

20.36 

0.65 


48 

304.2 ± 1.9 

Y 

4921 

1.73 

-1.89 ±0.13 

189086 

12 21 41.79 

-32 04 15.9 

21.46 

-0.21 


4 

315.5 ± 23.5 

Y 




194325 

12 21 42.51 

-31 59 10.0 

21.04 

0.51 


24 

301.9 ± 3.0 

Y 

5171 

2.15 

-2.76 ±0.43 

197616 

12 21 42.62 

-31 56 44.2 

21.11 

0.36 


15 

94.9 ± 6.6 

N 




201098 

12 21 42.78 

-31 54 26.7 

20.21 

0.82 


22 

-85.0 ± 2.3 

N 




196797 

12 21 42.79 

-31 57 04.2 

22.37 

0.43 


6 

300.1 ± 21.8 

Y 




195247 

12 21 44.45 

-31 58 21.0 

21.92 

0.48 


11 

317.5 ± 16.8 

Y 




196052 

12 21 44.57 

-31 57 37.6 

21.85 

0.56 


12 

303.7 ± 4.1 

Y 




200162 

12 21 44.65 

-31 55 17.7 

19.34 

0.98 


79 

68.1 ± 1.5 

N 




194736 

12 21 44.65 

-31 58 47.6 

21.95 

0.49 


10 

288.3 ± 14.5 

Y 




193869 

12 21 45.96 

-31 59 34.4 

20.11 

0.50 


46 

237.3 ± 2.1 

N 




183842 

12 21 46.10 

-32 07 02.8 

19.79 

0.67 


55 

-19.6 ± 1.7 

N 




202029 

12 21 47.09 

-31 53 36.0 

21.69 

0.50 


6 

10.1 ± 9.0 

N 




197129 

12 21 48.34 

-31 56 52.2 

21.50 

0.62 


17 

315.7 ± 3.7 

N 




192059 

12 21 48.39 

-32 01 21.8 

22.36 

0.59 


7 

172.4 ± 10.0 

N 





Pisces II 


9004 

22 58 17.52 

±05 55 17.5 

20.29 

0.73 


58 

-224.9 ± 1.6 

Y 

4787 

1.34 

-2.38 ±0.13 

9618 

22 58 20.35 

±05 58 08.8 

22.54 

0.59 


8 

-301.6 ± 6.5 

N 




9833 

22 58 21.22 

±05 57 20.3 

21.76 

0.61 


19 

-226.9 ± 3.2 

Y 




10215 

22 58 22.88 

±05 57 35.5 

21.53 

0.65 


21 

-25.9 ± 3.1 

N 




10694 

22 58 25.06 

±05 57 20.4 

19.94 

1.05 


59 

-232.0 ± 1.6 

Y 

4132 

0.81 

-2.70 ±0.11 

11592 

22 58 28.74 

±05 56 56.9 

19.91 

1.08 


82 

4.3 ± 1.5 

N 




12924 

22 58 34.10 

±05 57 43.6 

21.19 

0.63 


31 

-221.6 ± 2.9 

Y 

4955 

1.82 

-2.10 ±0.18 

13387 

22 58 35.91 

±05 57 22.4 

22.30 

0.49 


11 

-215.8 ± 7.6 

Y 




13560 

22 58 36.71 

±05 58 06.2 

21.48 

0.58 


23 

-232.6 ± 5.3 

Y 

5049 

1.99 

-2.15 ±0.28 

13757 

22 58 37.49 

±05 56 24.8 

21.25 

0.64 


27 

-102.1 ± 2.3 

N 




14179 

22 58 39.36 

±05 56 00.7 

22.42 

0.54 


8 

-224.8 ± 9.6 

Y 




16716 

22 58 51.19 

±05 59 57.1 

20.33 

0.74 


25 

-5.1 ± 10.9 

N 




17500 

22 58 54.78 

±05 57 47.5 

20.56 

0.83 


25 

83.7 ± 2.2 

N 





Laevens 1 


302 

11 36 13.55 

-10 53 02.9 

20.94 


1.29 

20 

149.8 ± 2.9 

Y 

4977 

1.53 

-1.59 ±0.15 

374 

11 36 13.76 

-10 52 48.5 

19.90 


1.32 

68 

151.9 ± 1.7 

Y 

4966 

1.09 

-1.66 ±0.11 

420 

11 36 13.90 

-10 52 27.3 

19.44 


1.59 

105 

149.8 ± 1.5 

Y 

4626 

0.69 

-1.55 ±0.11 

378 

11 36 14.53 

-10 52 43.7 

21.55 


1.22 

19 

153.5 ± 3.1 

Y 




93 

11 36 15.91 

-10 52 43.6 

19.68 


1.52 

88 

147.2 ± 1.6 

Y 

4726 

0.84 

-1.65 ±0.11 

1710 

11 36 16.52 

-10 52 46.7 

19.44 


1.31 

16 

143.6 ± 4.5 

Y 




1715 

11 36 16.59 

-10 52 46.2 

18.91 


0.85 

66 

72.0 ± 1.8 

N 




367 

11 36 16.89 

-10 52 43.7 

22.54 


1.14 

8 

153.8 ± 10.0 

Y 




1717 

11 36 17.26 

-10 52 46.2 

18.99 


0.70 

47 

266.0 ± 2.2 

N 




1972 

11 36 17.70 

-10 52 08.0 

22.17 


1.04 

10 

141.4 ± 4.7 

Y 




1997 

11 36 18.37 

-10 52 00.7 

22.60 


1.07 

6 

166.6 ± 11.5 

N 




1684 

11 36 18.59 

-10 52 49.7 

21.23 


1.23 

21 

150.9 ± 3.1 

Y 

5064 

1.70 

-2.10 ±0.21 

399 

11 36 19.58 

-10 52 37.1 

19.15 


0.70 

80 

155.3 ± 1.8 

N 




963 

11 36 20.11 

-10 52 38.2 

18.90 


1.87 

150 

149.2 ± 1.5 

Y 

4331 

0.23 

-1.78 ±0.11 


^ To convert to S/N per pixel, multiply by 0.57. 
^ Non-member based on CMD position. 
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Figure 4. The distributions of radial velocities for stars observed 
in the field of each satellite. Stars in the red portion of the his¬ 
tograms are deemed members. The dashed line shows (nheiio )5 
N is the number of member stars. The membership cut includes 
99% of member stars (2.58f7i, from the mean velocity). Stars whose 
velocity error bars overlap any part of the allowed velocity range 
are considered members. 

lated velocity dispersions as described in Section [5^ as¬ 
suming SsysV = 2.2 km s“^, the original value calculated 
by ISimon fc Gehd (|2007ll . instead of 1.49 km s“^. The 
increased systematic error would make the upper lim¬ 
its on the velocity dispersions about 5% less stringent, 
and it would reduce the measured velocity dispersion of 
Pisces II by 13%, well within the uncertainty. In sum¬ 
mary, our lower value of SsysV than previous works does 
not change our conclusions. 

Table [5] presents the velocities for stars in the satel¬ 
lite galaxies. The table excludes stars where velocity 
errors exceeded 30 km s“^ and stars where velocity mea¬ 
surement was not possible. Figure |4] shows velocity his¬ 
tograms for each satellite in bins of 2.5 km s“^. The 
domain of the plots is 100 km s“^ centered on the mean 
velocity for each satellite. In Section [51 we use these 
velocities to measure the kinematic properties of each 
satellite. 


4.2. Metallicity Measurements 

We measured metallicities by fitting synthetic spectra 
to the obser v ed spe c tra fo llowing the same procedure as 
iKirbv et al.1 (l2008al 1201011 . First, we normalized each 
spectrum to its continuum. To do so, we divided the 
spectrum by a spline fit to spectral regions generally free 
of absorption lines. Then, we matched the spectrum 
to a grid of synthe tic spectra ge nerated with ATLAS9 
model at mospheres (lKuruczlll993[l and the synthesis code 
MOOG ([Snedenl 119731 1. IW measured iron abundances, 
[Fe/H]3, by fitting only to iron absorption lines and ig¬ 
noring other spectral regions. Please refer to IKirbv et al.l 
(1201011 for more detail. 

We estimated the random uncertainty on [Fe/H] from 
the covariance matrix of the fit to the synthetic spec¬ 
tral grid. There is also a systematic uncertainty of 
0.11 dex (IKirbv et al.ll2010ll . The final error bar is the 
quadrature sum of the random and systematic uncertain¬ 
ties. 

Table [2] gives effective temperatures (Tes), surface 
gravities (log 5 ), and [Fe/H] measurements for member 
stars with (5[Fe/H] < 0.5 and S/N > 20 A“^. Fig- 

® We adopted the solar abundances of lAndersfcGreyessd 1 1989|) 
except for iron: 12 + log(n(Fe)/n(H)) = 7.52 IKjneden et al.lll892li . 


Figure 5. The distribution of [Fe/H] in each satellite. There 
are fewer stars in this figure than in Figure [^because we could not 
measure [Fe/H] for all radial velocity members. N is the number of 
stars for which [Fe/H] measurements with errors less than 0.5 dex 
were possible. 
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Figure 6. DEIMOS spectra of Lae vens 1 stars 420 (black) and 93 
(red), which were also observed bv (Bonifacio et al.l II2015I ). The 
bold (black and red) spectra show the spectra as we observed 
them. The faded (gray and pink) spectra show the spectra as 
they would hav e appeared at the radial velocities measured by 
IBonifacio et al.l The visible separation between the faded spectra 
shows that our spectra have the precision to rule out a velocity 
separation of 10.2 km s~^. 

ure [5] shows the distribution of these measurements in 
bins of 0.25 dex. By nature, ultra-faint satellites have 
only handfuls of red giants bright enough for metallic¬ 
ity measurements. As a result, the metallicity distribu¬ 
tions in Figure [5] have only 4-6 stars. Accurate quan¬ 
tification of metallicity distributions and chemical evo¬ 
lution requires larger samples. Nonetheless, we present 
rudimentary metallicity averages and dispersions in Sec¬ 
tion jB] 


4.3. Comparison to Other Measurements 

IBonifacio et al.l (I2015D obtained VLT/X-Shooter spec¬ 
tra of two stars in Laevens I. Both stars are also 
in our spectroscopic sa mple. We call them 420 and 
93, and IBonifacio etldl call them JII36I3—105227 and 
J113615—105244. They measured a velocity difference 
between the two stars of 10.2 ± 5.7 km s“^, whereas we 
measured a difference of only 2.6 ± 2.2 km s“^. 

Figure [ 6 ] shows our DEIMOS spectra in the observer 
frame. The bold black and red spectra show the spec¬ 
tra as we observed them, and the faded (gray and pink) 
lines show the spectra as t hey would have appeared at 
the velocities measured by IBonifacio et al.l (|2015D . The 
bold spectra are nearly superposed on top of each other, 
whereas the faded spectra are visibly separated. Thus, 
the figure shows that our spectra have the precision to 
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discern between a 10.2 km and a 2.6 km differ¬ 
ence. 

There are several possible reasons for the discrepancy 
between our two studies. First, one or both of the stars 
could have a variable radial velocity. In that case, we 
happened to observe the two stars when they had nearly 
the same radial velocity, and our inability to resolve the 
velocity dispersion (Section 15.21) is unaffected. Second, 
one of our studies could be subject to a systematic error. 
For example, we corrected for slit mis-centering by ap¬ 
plying a velocity shift based on telluric absorption. (The 
wavelengt h scale in Figur e [5] includes the telluric cor¬ 
rection.) iBonifacio et al.l were unable to apply a sim¬ 
ilar correction. Therefore, one possible systematic er¬ 
ror in their study was slightly different positions of the 
two stars perpendicu lar to the X-Shooter slit. Regard¬ 
less, IBonifacio et al.l state that the velocity dispersion of 
Laevens 1 is consistent with zero when the measurement 
errors are taken into account. This basic result agrees 
with our own conclusions ISection 15.21) . 

We measured [Fe/H] = —1.55 ± 0.1 1 and —1.65 ± 
0.11 for stars 420 and 93, respectively. IBonifacio et al.l 
measured [Fe/H] = —1.73 ± 0.26 and —1.67 ± 0.28 
from an equivalent w idth analysis using MyGIsFOS 
(|Sbordone et al.ll2014fl . Our measurements are entirely 
consistent within the measurement uncertainty. 

5. KINEMATICS 
5.1. Membership 

Ultra-faint satellites have very low surface brightness. 
In some cases, most of the point sources in the vicinity 
of an ultra-faint satellite—even at the center—are fore¬ 
ground dwarf stars. We avoided most of the contamina¬ 
tion by selecting stars with colors and magnitudes ap¬ 
propriate for red giants at the distance of each satellite 
(Section [2] and Figure [T]). We removed the remainder of 
the contaminants by imposing a radial velocity cut. 

We determined membership in each satellite by using 
its mean velocity, (uheiio), and velocity dispersion, Uy. 
First, we started with guesses for (uheiio) and tr^. We re¬ 
quired member stars to have juheiio — (fheiio)| < 2.58tT«. 
For a Gaussian velocity distribution, this cut includes 
99% of member stars. Because for most of the satel¬ 
lites is very small, a hard cut would exclude low-S/N 
stars whose measured velocities are discrepant by more 
than 2.58(71, simply because their uncertainties are larger 
than that. Therefore, we extend the membership crite¬ 
rion to any star whose Icr error bar overlaps the velocity 
range for member stars. We determined (uheiio) and 
from the list of member stars following the procedure in 
Section O Then, we used these new values to reform 
the member list. We repeated this procedure until the 
member list did not change from one iteration to the 
next. 

Foreground dwarfs can also be identified by spectral 
features sensiti ve to surface gravity, such as the Na i dou¬ 
blet at 8190 A (iS oinrad fc Tavloil[l971 :lCo henlll978D and 
Mg I 8807 (jBattaglia fc Starkenburd 2012D . We found 
several stars with very strong Na doublets and Mg i lines 
among the Hydra H and Pisces H samples, but they were 
already excluded by the radial velocity cut. Hence, we 
did not need to impose any additional membership crite¬ 
ria. We discuss Mg i again in Section [5i2] in reference to a 


star that barely missed the membership cut in Hydra H. 

Two of the three blue stars that we targeted in 
Laevens 1 are non-members on the b asis of their ra¬ 
dial velocities. iBelokurov et al.l (|2014fl speculated that 
these stars might be blue loop stars. If they were, then 
Laevens 1 must have a young stellar population, which 
would argue strongly that it is a star-forming galaxy, 
not a GC. However, the non-membership of these two 
stars negates most of the evidence that Laevens 1 has a 
very y oung stellar population, although IBonifacio et al.l 
(|2015l) argued that the blue extension of the main se¬ 
quence might indicate the presence of a ~2 Gyr old 
population. The third blue star, 399, has Uheiio = 
155.3±1.8 km s“^, which is 6.0 km s“^ from the mean ra¬ 
dial velocity. This star formally passes the velocity mem¬ 
bership cut. However, in addition to its unusual color, 
it has a more discrepant radial velocity than any of the 
other member stars, and it is farther from the center of 
the cluster than all but one confirmed member. We con¬ 
sidered the possibility that this star is a Cepheid mem- 
ber of Laevens 1. However, the Pan-STARRSl survey 
(|Kaiser et al.ll2010ll obtained six observations of this star 
in each of the rpi, zpi, and ypi bands, and the rms scat¬ 
ter of these observations is consistent with measurement 
noise (^ 0.03 mag, B. Sesar, priv. comm.). Since Popula¬ 
tion H Cepheids have light curve amp litudes between 0.5 
and 1 .2 mag in the Johnson i?—band (jWallerstein fc Co^ 
[TMl . we conclude that the blue star is not a Cepheid 
variable. We discuss the impact that including star 399 
as a member would have on Uy in Section o 

5.2. Mean Velocities, Velocity Dispersions, and Masses 

We estimated (uheiio) and g„ for the thre e satellites 
in the same manner that iKirbv et al.l (l2013all measured 
these values for th e UFD Segue 2 . Th at method, in 
turn, was based on iWalker et all ’s (|2006( ) procedure for 
measuring the velocity dispersion of the Fornax dwarf 
spheroidal galaxy (dSph). The method uses maximum 
likelihood statistics and a Monte Carlo Markov chain 
(MCMC). We maximized the logarithm of the likelihood 
(L) that the given values of (uheiio) and Uy described the 
observed velocity distribution, including the uncertainty 
estimates for individual stars. 


logL = ^ ^ (log(((5uheiio)? + (tI) 

i 

I 1 f (('ahelio)i ~ ('Chelio))^ \ , 

^2^^, (<5uhelio)?+a2 ) 

As initial guesses for (uheiio) and (t„, we started with 
the mean and standard deviation of Uheiio over the ve¬ 
locity range shown in Figure ID The final results are not 
sensitive to these guesses. We explored the parameter 
space with a Metropolis-Hastings implementation of an 
MCMC. The length of the chain was 10^ trials. 

As a test of our procedure, we measured (uheiio) = 
— 144.4 ± 1.3 km s“^ and cr„ = q km s“^ for 

M22 from the 2009 Oct 1 4 obse rvation. In compari¬ 
son, iPeterson CudworthI (|1994[ ) measured (uheiio) = 
—148.8 ± 0.8 km s“^ and a ,, = 6.6 ± 0 .8 km s“^ within 
7' of the cluster center, and lLane et al.l (|2009ll measured 
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Figure 7. The probability distributions for The histograms 
show successful MCMC trials. Vertical lines show the measured 
velocity dispersion (Pisces II) or upper limits (Hydra II and 
Laevens 1). 



(wheiio) = —144.9 ± 0.3 km s ^ and a central dispersion 
of Oy = 6.8 ± 0.9 km s ~^. Our me asurement of (iiheiio) is 
consistent with that of iLane et al.l. and our measurem ent 
of (Jy is con sistent with both IPeterson fc CudworthI and 
ILane et ahl The velocity dispersion in M22 decreases to 
abou t 5 km s“^ at 8' from the cluster center (|Lane et al.1 
I2009H . However, our data is largely insensitive to the de¬ 
cline is Gy because the maximum radial extent of our 
spectroscopy is 8', and most of our M22 targets are 
within 5'. 

Table 13] gives (uheiio) and cjy with errors. It also gives 
velocities relative to the Galactic standard of rest (GSR) 
assuming that the Sun’s orbital velocity is 220 km s“^. 
Hydra H is receding from the Galactic center. There¬ 
fore, it is past its pericenter on the way to its apocen- 
ter. Its heliocentric radial velocity is also similar to that 
of the gas in the leading arm of the Ma g ellanic stream 
(1 Put man et al.1 II998I : iBriins et al.l 120051 : iNidever et al.l 
1200811 . This finding strengthens the po tential for Hy¬ 
dra H to be associated with the LMC (IKoposov et al.l 
l2015aHBechtol et al.ll2015l: iDeason et al.ll2015ll . Pisces H 
is approaching the Galactic center, which means that it 
is on the way to its pericenter. Finally, Laevens I has a 
very small uqsR) meaning that it is at pericenter or apoc- 
enter. Because it is so distant (at least 140 kpc from the 
Galactic center), it is either close to apocenter or on a 
quasi-circular orbit. 

Figure |7| shows the distribution of ay for accepted 
MCMC trials. These distributions are equivalent to 
the probability distribution for cr^,. The distribution 
for Pisces H has a well-defined peak separated from 
zero. The dashed line shows the median value of ay = 
5.4 km s“^. We estimated asymmetric error bars by cal¬ 
culating the values of ay that bracket 68.3% of the prob¬ 
ability distribution. 

The distributions for Hydra H and Laevens I are con¬ 
centrated near zero. Hence, our measurements cannot 
resolve ay for these two satellites. We estimated two sets 
of upper limits by finding the value of cr„ that exceeds 
90% and 95% of the MCMC trials. Table |3I gives both 
upper limits for both satellites. 

Hydra H and Laevens I have three stars with some¬ 
what ambiguous membership. Stars I97I29 and 188838 
in Hydra H have velocities larger than (wheiio) by 12.6 
and 15.9 km s“^. If our membership cut were more in¬ 
clusive than 2.58ay, then these stars could be considered 


members. Unfortunately, neither star has a visible Na i 
doublet to help define its status. The Mg i equivalent 
width (EW) of star 188838 is 0.34 ±0.07 A, and the com¬ 
bined EW of of the two redder lines of the Ca ii infrared 
triplet (CaT) is 4.1 ± 0.3 A. These measurements fall 
right on the dividing line between dwarf s and giants de¬ 
fined by IBattaglia fc Starkenburd (I2012f) . Because most 
metal-poor dSph stars fall well below the dividing line, 
this test disfavors membership for star 188838. The Mg i 
line in star 197129 is too noisy to be useful. 

Including star 192179 as a member does not change 
our qualitative conclusions. Instead, the 90% C.L. upper 
limit on ay rises from 3.6 km s“^ to 6.0 km s“^. On 
the other hand, including star 188838 resolves the veloc¬ 
ity dispersion as ay = 6.4l^ g km s“^. Including both 
stars results in ay = 6.7^^ g km s“^. However, both 
stars would be > 2cr„ outliers even with the larger ay. 
Eurthermore, Figure [T] shows that both stars are redder 
relative to the M92 isochrone than any other member. 
Figure m also shows that star 188838 would be the most 
distant member of Hydra 11. Because the majority of 
the evidence disfavors membership, we ruled both stars 
as non-members. Including or excluding the two possi¬ 
ble HB stars, 189086 and 191385, changes the limits on 
ay by less than 2%. Based on their CMD positions, we 
ruled 189086 a member and 191385 a non-member. 

Including star 399 in Laevens 1 resolves the velocity 
dispersion as ay = 2.71^;® km s“^. As discussed in Sec- 
tion l5.ll the star is bright and blue. As a result, star 399 
is probably not a member. 

The total mass of a spherical stellar system in dynami¬ 
cal equil ibrium is related t o the square of the velocity dis¬ 
persion. IWolf et (|2010ll showed that the total mass is 
poorly constrained because any possible underlying dark 
matter has an unknown mass profile, and the velocity 
anisotropy cannot be determined from a small sample of 
radial velocities. However, the mass within the half-light 
radius, Mif 2 , is well-constrained: M 1/2 = 4G“^cr^r/i. 
Although Th is the 2-D half-light radius, the formula in¬ 
fers the mass enclosed within the 3-D half-light radius. 

Table HI gives log Mi /2 for Pisces H and upper limits 
for Hydra H and Laevens I, which have only upper lim¬ 
its for ay. Pisces H has a dynamic al mass on par with 
dwarf galaxies of similar luminosity (|Strigari et al.l[2008l : 
IWolf et ani20I0[) . For example, its lumi nosity and mass 
are v ery similar to Ganes Venatici H (|Simon fc Gehal 
l2007f) . The mass limits for Hydra H and Laevens 1 are 
not stringent enough to make these satellites unusually 
light. The 90% G.L. for Hydra H is M 1/2 < 1.0x 10® Mq, 
which is less than a factor of 2 smaller than Pisces H. The 
90% G.L. mass limit for Segue 2, the least massive galaxy, 
is 7 times more stringent than Hydra H. The limit for 
Laevens 1 is slightly more stringent than Hydra H, but 
w e argue below that Laeven s 1 is not a galaxy. 

iWillman fc Stradeil (1201211 proposed that a “galaxy” 
be defined as a stellar system that cannot be explained by 
a combination of baryons and Newton’s laws of gravity. 
The mass-to-light ratio for Pisces H is 370^*^240 -^©/ 
This value is far too large to be explained by baryons 
alone, even for the oldest stellar populations in the Uni¬ 
verse. In fact, 99% of the successful MGMG trials have 
M/L > 10 Mq/Lq. Hence, Pisces H satisfies the defi¬ 
nition of a galaxy. However, upper limits on the mass- 
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Table 3 

Satellite Properties 


Property 

Hydra H 

Pisces II 

Laevens 1 

-^member 

13 

7 

10 

log(Lv/LQ) 

3.90 ±0.10 

3.93 ±0.20 

3.65 ±0.08 

Vh (arcmin) 

-*-•^-0.2 

1.1 ±0.1 

0 47+0-04 

rh (pc) 

66lg^ 

58 ±7 

19 ±2 

(Vhelio) (km s-l) 

303.1 ± 1.4 

-226.5 ±2.7 

149.3 ± 1.2 

■CGSR (km s-l) 

135.4 

-79.9 

4.6 

(Tv (km s“^) 

< 3.6 (90% C.L.) 

5.4t3-6 

< 3.9 (90% C.L.) 


< 4.5 (95% C.L.) 

6.2/°;3 

< 4.8 (95% C.L.) 

log(iW'i/2/AL0) 

< 5.9 (90% C.L.) 

< 5.5 (90% C.L.) 

< 6.1 (95% C.L.) 

'^'*^-240 

< 5.6 (95% C.L.) 

{Mq/Lq) 

< 200 (90% C.L.) 

< 130 (90% C.L.) 

< 315 (95% C.L.) 


< 192 (95% C.L.) 

([Fe/H]> 

-2.02 ±0.08 

- 2.45 ±0.07 

-1.68 ±0.05 

<r([Fe/H]) 

0 

0.48/°;™ 

< 0.40 (90% C.L.) 

< 0.53 (95% C.L.) 


Re f erences. — The measu reme nts of logLy and rn come from IMartin et al.l 
1 I 2 OI 5 I I , [Belokurov et aO II2010I 1 , and ILaevens et al.l II2014I 1 . 


^ Mass-to-light ratio within the half-light radius. 


to-light ratios on Hydra II and Laevens 1 do not help in 
deciding if they are galaxies. In Section El we use chemi¬ 
cal evidence to resolve the nature of these two satellites. 

One of the reasons dwarf galaxies are interesting is that 
they are targets for the detection of gamma rays due 
to dark matter self-annihilation. Pisces II has similar 
structural properties {rh and cr„) and distance to Canes 
Venatici II. Hence, the two galaxies’ potential for the 
detection of the gamma ray signal is about the same. 
iBonnivard et al.l (1201511 found that Canes Venatici H is 
not the most promising dwarf galaxy to search for self- 
annihilation, but it would contribute significantly to an 
analysis that stacks the Fermi g amma ray telescope ob¬ 
serva tions of all of the dwarfs (jAckermann et aLri20l4 
EqU, e.g.,). 

6 . METALLICITY 





0.0 0.5 1.0 1.5 2.00.0 0.5 1.0 1.5 2.00.0 0.5 1.0 1.5 2.0 
o([Fe/H]) C5([Fe/H]) a([Fe/H]) 

Figure 8. The probability distributions for (T([Fe/H]). The his¬ 
tograms show successful MCMC trials. Vertical lines show the 
measured metallicity dispersions (Hydra II and Pisces II) or upper 
limits (Laevens 1). 


The stellar mass of a UFD is insufficient to retain su¬ 
pernova ejecta. Nonetheless, all UFDs studied in suf- 
hcien t detail show evidence for metallicity dispersions 
(e.g., Frebel et all 12(3101: iNorris et 311120101 : iKirbv et al.l 
l2013anVargas et al.ll2013lb Their ability to self-enrich 
with iron implies that they have or once had enough 
mass to prevent met al-enriched supernova ejec ta from 
escaping. As a result. IWillman fc Stradeil (j2012D consid¬ 
ered a dispersion in [Fe/H] sufficient for classification as 
a galaxy for most stellar systems 0 
We calculated the metallicity mean and dispersion in 
the same way that we computed (uheiio) and a^. In anal¬ 
ogy to Equation [1] we maximized the likelihood that the 
metallicity distribution has a mean ([Fe/H]) and disper¬ 
sion cr([Fe/H]) as follows: 

log£= ^log(2^) + i f] (log(<5[Fe/H])? + a([Fe/H])^) 

i 

^ Some massive GCs, like to Centauri JNorris Sz Da Costalll995l^ . 
have a dispersion in metallicity but no kinematic evidence for dark 
matter. 


1 ^ f ([Fe/H], - ([Fe/H]))^ \ 

+ 2^G<5[Fe/H])|+a([Fe/H])2; ^ ^ 

We used the metallicity measurements shown in Table O 
i.e., those with (5[Fe/H] < 0.5 and S/N > 20 A“^. The 
values were determined through lO"^ MCMC trials. 

Figure El shows the distribution of (T([Fe/H]) for the 
successful MCMC trials. Whereas the distributions for 
Hydra H and Pisces H are somewhat separated from zero, 
the distribution for Laevens 1 piles up at zero. Hence, 
we marginally resolve the metallicity dispersion for Hy¬ 
dra H and Pisces H: 79% and 85% of the MCMC trials 
have crpe/H] > 0-2, respectively. On the other hand, we 
measured only an upper limit for Laevens 1. These val¬ 
ues are shown in Table El 

Although we did not resolve the velocity dispersion for 
Hydra H, we did detect a metallicity dispersion, albeit 
with somewhat low significance. Because it seems that 
Hydra H enriched itself with iron, we tentatively assign it 
a classification as a galaxy. This spectroscopic classifica¬ 
tion supports the phot ometric classificati on based on its 
large half-light radius ([Martin et al.l[2M^ . Pisces H also 
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Figure 9. The luminosity—metallicity relation for Milky Way 
satellite galaxies. Hydra II and Pisces II fall in the midst of other 
UFDs. Pisces II falls almost directly on top of Leo IV. Like other 
GCs, Laevens 1 does not conform to the relation. The black points 
are ([Fe /H]) measurements for Milky Way dSphs from IKirbv et ahl 
Il2013bl ). The dashed line is a least-squares fit to the black points 
excluding Segue 2, the least luminous galaxy on the plot. The 
dotted lines show the rms dispersion about the relation. 

has a marginally resolved metallicity dispersion, consis¬ 
tent with its dynamical classification as a galaxy. 

Dwarf galaxies obey a tight rel ationship between lu¬ 
minosity and average metallicity (iSkillman et aH Il989t 
iMateol 119981 lKirbvll20i^ iKirbv et al.112013^ Hi fact, 
the luminosity-metallicity relation (LZR) can be used as 
a diagnostic for whether a stellar system is a ga laxy. Fig¬ 
ure ini shows the LZR from IKirbv et al.l (|2013bf) with the 
three new satellites. Hydra H and Pisces H lie in the 
same region of the diagram as other UFDs, like Leo IV 
and Canes Venatici H. Thus, the LZR supports their 
identification as galaxies. 

On the other hand, Laevens 1 is too metal-rich for its 
luminosity (or too faint for its metallicity). GCs do not 
obey any LZR. Tidal dwarf galaxies also do not obey the 
LZR, but tidal d warfs tend to be close to solar metallicity 
l|Duc et al.ll2000H . Therefore, Laevens I’s position in Fig- 
ure[9]suggests that it is a GC or a severely tidally stripped 
dwarf galaxy, in which the removal of stars has reduced 
Laeven s I’s mass but not its metallicity. IKirbv et al.1 
(|20I3all proposed this scenario for Segue 2, but Segue 2 
has a metallicity dispersion. We detected no metallicity 
dispersion in Laevens I, favoring its identification as a 
GC. 

Laevens 1 has four pieces of evidence that suggest that 
it is a GC, not a galaxy. F irst, it has a smaller half-light 
radius than known UFDs (|Laevens et al.ll^Oldf) . Second, 
it does not have strong evidence for dark matter {M/L < 
130 Mq/Lq, 90% C.L.). Third, its metallicity dispersion 
is less than 0.40 dex (90% C.L.). Fourth, it does not 
obey the LZR for dwarf galaxies. Additionally, most of 
the blue stars that may have suggested that Laevens 1 
has a young stellar population turned out to be probable 
non-members ISection IS.ip . The various lines of evidence 
taken together favor a cluster classification rather than 
a dwarf galaxy. 


Star 10694, the brightest member of Pisces H, is 
carbon-rich. Figure |3] shows CN absorption between 
8300 A and 8400 A. Bright red giants should be de¬ 
stroying carbon, not creating it. Star 10694 could be an 
AGB star that is currently dredging up carbon. Alter¬ 
natively, it could have acquired carbon from a recently 
defunct AGB companion. Its large luminosity favors the 
scenario that it is itself an AGB star. 


7. SUMMARY 

We obtained Keck/DEIMOS spectroscopy for three 
newly discovered Milky Way satellites: Hydra H, 
Pisces H, and Laevens 1. We measured radial veloci¬ 
ties and metallicities for individual candidate member 
stars. We identified 13, 7, and 10 member stars in the 
three satellites, respectively, on the basis of their radial 
velocities. Most of the members are red giants, although 
there is one HB and one AGB star in Hydra H and one 
carbon-rich star in Pisces H. 

We could not resolve the velocity dispersion of Hy¬ 
dra H, but we did measure a non-zero dispersion in metal¬ 
licity (cr([Fe/H]) = 0.40to'26)' Because it seems to have 
chemically enriched itself in iron and because it has a 
large half-light radius (66 pc, iBelokurov et al.l[2010ll , Hy¬ 
dra H is more likely to be a dwarf galaxy than a GG. It 
also has a radial velocity similar to the gas in the leading 
arm of the Magellanic stream at its location. Therefore, 
it may have fallen into the Milky Way with the Magel¬ 
lanic Glouds. 

Pisces H is a bona fide galaxy that inhabits a massive 
dark matter halo. We measured a velocity dispersion in 
Pisces H that is far in excess of what would be expected 
on the basis of its stellar mass alone. We also marginally 
resolved a metallicity dispersion ((T([Fe/H]) = 0.48lo;2g), 
indicating that Pisces H is not only massive now, but it 
was also massive enough during star formation to retain 
supernova ejecta. 

Laevens 1 is more likely a GG than a galaxy. We did 
not resolve a dispersion in velocity or metallicity, and 
Laevens 1 does not obey the LZR for dwarf galaxies, 
whereas Hydra H and Pisces H do. Although the upper 
limit on the mass-to-light ratio {M/L < 130 Mq/Lq, 
90% G.L.) does not rule out that Laevens 1 could be dark 
matter-dominated, the low metallicity dispersion and the 
inconsistency with the LZR suggest that the system is a 
gl obular cluster. 

IBelokurov et al.l (|2014ll found several bright, blue stars 
in the vicinity of Laevens 1. They suggested that these 
could be blue loop stars, which would signify the presence 
of a very young stellar population in Laevens 1. However, 
we found these stars to be probable non-members on the 
basis of radial velocity. 

With this work, we have spectroscopically confirmed 
two dwarf galaxies and one GG. The dwarf galaxies are 
especially interesting for their application to dark mat¬ 
ter physics. In particular, Pisces H has a mass-to-light 
ratio within the half-light radius of 370;1^24o 
though this does not make it one of the most dark matter- 
dominated galaxies known, it is still massive enough 
to warrant inclusion in a search for gamma rays due 
to dark matter self-annihilation ( e.g., lAckermann et all 
IMlIMlfBonnivard et al.lIM^ . 
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